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In a recent Letter by the authors [I.S. Burmistrov and A.M.M. Pruisken, Phys. Rev. Lett. 101, 
056801 (2008)] it was shown that single-electron devices (single electron transistor or SET) dis- 
play "macroscopic charge quantization" which is completely analogous to the quantum Hall effect 
observed on very different electron systems. In this investigation we present more detail on these 
new findings. Based on the Ambegaokar-Eckern-Schon (AES) theory of the Coulomb blockade we 
introduce a general response theory that probes the sensitivity of SET to changes in the boundary 
conditions. This response theory defines a new set of physical observables and we establish the 
contact with the standard results obtained from ordinary linear response theory. The response pa- 
rameters generally define the renormalization behavior of the SET in the entire regime from weak 
coupling with large values of the tunneling conductance all the way down to the strong coupling 
phase where the system displays the Coulomb blockade. We introduce a general criterion for charge 
quantization that is analogous to the Thouless criterion for Anderson localization. We present the 
results of detailed computations on the weak coupling side of the theory, i.e. both perturbative and 
non-perturbative (instantons). Based on an effective theory in terms of quantum spins we study the 
quantum critical behavior of the AES model on the strong coupling side. Consequently, a unifying 
scaling diagram of the SET is obtained. This diagram displays all the super universal topological 
features of the 9 angle concept that previously arose in the theory of the quantum Hall effect. 

PACS numbers: 73.23.Hk, 73.43.-f, 73.43.Nq 



I. INTRODUCTION 
A. Coulomb blockade 

The Coulomb blockade in nanostructures is one of the 
cornerstones of modern condensed matter physics. The 
simplest approach to electron tunneling through quan- 
tum dots was proposed by Ambcgaokar, Eckern and 
Schon in 1982J-2 Their model (in brief AES model) be- 
came the focus of a stream of experimental and theoreti- 
cal papers^^ following the first experimental indications 
of "macroscopic charge quantization" in single electron 
devices in 1991 

To experimentally control the transport of electrons^ 
one generally uses the so-called "single electron transis- 
tor" or SET4 This is a mesoscopic metallic island that 
is capacitively coupled to a gate and connected to two 
metallic reservoirs through tunneling contacts with a to- 
tal conductance g (see Fig. QJ,). 

The experimental conditions of the AES model are lim- 
ited and extremely well known^ ' 10 ! 11 The model is nev- 
ertheless richly complex and much of the physical conse- 
quences have remained unknown. Over the years, how- 
ever, it has slowly become more evident that the AES 
theory of the Coulomb blockade is in many ways similar 
to the theory of the quantum Hall effect^ For example, 
the AES model is asymptotically free in + 1 space-time 
dimension, possesses instantons and has an instanton an- 
gle 9. This immediately raises the fundamental question 
whether the experimental phenomenon of "macroscopic 
charge quantization" in the SET is possibly related to the 



"robust quantization" of the Hall conductance observed 
on very different electronic systems. 

The AES model has a number of very significant ad- 
vantages as compared to the more conventional theories 
of the 9 vacuum or instanton vacuum. For example, the 
winding numbers of the theory ( "topological charge" ) are 
quantized at the outset of the problem. This is quite un- 
like the usual situation where the historical controversies 
in quantum field theory continue to haunt the subject. 
For example, it has been pointed out only very recently 
that the 9 vacuum concept generally displays "massless 
chiral edge excitations" that are very different from those 
in the "bulk" of the system. Disentwining these different 
types of excitation is synonymous for separating the frac- 
tional topological sectors of the theory from the integral 
ones. 

Remarkably, it turns out the existence of "massless chi- 
ral edge excitations" in the problem automatically reveals 
the existence of the quantum Hall effect. This fundamen- 
tal phenomenon previously remained concealed. How- 
ever, it provides the resolution to longstanding problems 
such as the quantization of topological charge, the mean- 
ing of instantons and instanton gases etc. etc. 

The existence of "massless chiral edge excitations" has 
furthermore led to the idea of "super universality" which 
states that all the fundamental features of the instanton 
vacuum concept are precisely those of the quantum Hall 
effect^ These include not only the robust quantization 
of the Hall conductance but also the existence of "gapless 
excitations" at 9 = tt or, equivalently, "quantum critical- 
ity" of the quantum Hall plateau transition. 

R is of interest to know whether these new advances 



2 



possibly also apply to the AES model. In this case, the 
microscopic origins of the integral and fractional topo- 
logical sectors are far more obvious. For example, the 
integral sectors directly emerge from quantum statistics 
and they describe the quantum system (SET) in thermal 
equilibrium. On the other hand, the fractional topologi- 
cal sectors do not describe "edge" excitations but, rather, 
they have the meaning of perturbing external fields that 
take the SET out of thermal equilibrium. The great ad- 
vantage of the AES model, however, is that the 9 depen- 
dence can be studied on the strong coupling side. The 
AES model is therefore an outstanding laboratory where 
the various different aspects of "super universality" can 
be explored and investigated in great detail. 

It should be mentioned that the AES model in a differ- 
ent context is also known as the "circular brane model."— 
It is furthermore of direct physical interest in the theory 
of granular metals at intermediate temperatures^ 



B. Charge quantization 

The phrase "macroscopic charge quantization" usually 
refers to the charge of an isolated island that is discon- 
nected from the reservoirs. It is given by the naive strong 
coupling limit of the AES model where the tunneling con- 
ductance g is put equal to zero. 

This naive approach leads to the electrostatic picture 
of the Coulomb blockade where the average charge (Q) 
on the island is robustly quantized in units of e as the 
temperature (T) goes to absolute zero. This quantization 
breaks down for very special values of the gate voltage 

VW = e(k+l/2)/C g (1) 

where fc is an integer and C g denotes the gate capaci- 
tance. At these very special values of V g a first order 
quantum phase transition occurs separating two different 
phases with Q = k and k + 1 respectively. 

The electrostatic picture of the Coulomb blockade gets 
fundamentally complicated when the tunneling conduc- 
tance g is finite. It is well known, for example, that due 
to the strong charge fluctuations in the SET the aver- 
aged charge Q on the island is generally wn-quantized.— 
Despite the impressive list of existing theoretical work on 
both the strong coupling sido^^ (j<1) and weak cou- 
pling sid&i^ (g 3> 1) of the problem it is not known what 
the electrostatic or semi classical picture of the SET ex- 
actly stands for. This fundamental drawback clearly up- 
sets the concept of "robust charge quantization" in single 
electron devices. 



C. Outline of this investigation 

The main objective of this investigation is to show 
that the SET displays macroscopic charge quantization 
in much the same way as the two dimensional electron 



gas displays the quantum Hall effect. We develop a com- 
plete quantum theory of the SET and introduce a uni- 
fying scaling diagram that spans the entire range from 
weak to strong coupling. 

We benefit from the advances made over the years in 
the theory of the quantum Hall effect. We introduce, in 
particular, the "physical observables" of the AES theory 
that measure the sensitivity of the SET to changes in the 
boundary conditions. In the present context this means 
that the quantum system is taken out of thermal equi- 
librium by perturbing fields. The main problem to be 
solved is how to lay the bridge between the sensitivity to 
the boundary conditions on the one hand, and the stan- 
dard expressions for linear response obtained from the 
Kubo formalism on the other. 



1. Electrostatic picture revisited 

To start we briefly review the microscopic origins of 
the AES model and summarize the results known from 
previous work in Section HT1 

To see the concept of "physical observables" at work we 
consider in Section|TTT]the trivial case of an isolated island 
at finite T obtained by putting the tunneling conductance 
g equal to zero. This simple but instructive example sets 
the stage for most of the analysis in the remainder of this 
paper. 

We point out, first of all, that the averaged charge Q 
on the isolated island is a measure of the sensitivity of 
the system to changes in the boundary conditions. This 
notion immediately suggests a generalized Thouless cri- 
terion that relates the robust quantization of Q on the 
island to the appearance of an energy gap. 

Secondly, we show how the rcnormalization behavior of 
Q at finite T provides a complete knowledge of the low 
energy dynamics of the isolated island. This behavior 
involves two different kinds of fixed points, i.e. stable 
ones at Q = k which describe the robust quantization of 
charge as T goes to zero, and unstable ones at Q = fc+1/2 
describing the transition between the states Q = k and 
k + 1 of the island. 



2. Two sets of physical observables 

Armed with the insights obtained from the isolated is- 
land we next embark on the general problem with finite 
g in Section HVl We introduce two slightly different but 
physically equivalent sets of response parameters g' and 
q' . The different expressions that we obtain stem from 
slightly different ways of handling the fractional topolog- 
ical sectors of the AES theory. 

The first set is the simpler one which is a direct gen- 
eralization of the results obtained for an isolated island. 
The second set is slightly more involved but our findings 
permit a direct comparison with the expressions obtained 
from linear response theory. 
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In both cases, however, one may think of g' and q' 
in terms of the sensitivity of the SET to changes in the 
boundary conditions. In both cases also one may think 
of g' in terms of the SET conductance. The quantity q' is 
new and in general very different from the conventionally 
studied averaged charge Q on the island. Within linear 
response theory we express q' in terms of the antisym- 
metric current-current correlation function. Wc identify 
this new quantity with the previously unrecognized quasi 
particle charge of the SET. 

In complete analogy with the theory of the quantum 
Hall effect, we relate the conditions for "macroscopic 
charge quantization" g' = and q' = k with integer 
k to the appearance of an energy gap in the SET. For 
g = these conditions are identically the same as those 
obtained from the electrostatic picture of the SET. For 
finite g, however, these conditions describe an entirely 
different physical state of the SET. They describe the 
macroscopic quantization of the quasi particle charge, 
rather than the averaged charge on the island. 



3. Explicit computations 

This takes us to the second part of this investigation 
where we explicitly compute, in Sections |V] and IVI| 
the observable theory g' and q' in the various different 
regimes in g of interest. We benefit from having two 
different definitions of g' and q' . The different compu- 
tational schemes provide a direct check on the universal 
and non-universal parts of the AES theory. 

In Section [V] we consider the weak coupling phase g — > 
oo of the AES model. We report the detailed results for 
the renormalization group (3 functions based on ordinary 
perturbation theory as well as instantons. Even though 
this Section is self-contained, we refer the reader to the 
literature for a more detailed exposure to the instanton 
calculational technique^ 

In Section [VI] we address the strong coupling phase of 
the AES theory and study, in particular, the quantum 
critical behavior of the SET at finite g. For this pur- 
pose we first map the critical behavior of the AES model 
onto an effective theory of quantum spins. Wc employ 
Abrikosov's pseudo fermion technique and extract the (5 
functions of the AES theory near the critical point. 

The most important results of this investigation are en- 
capsulated in the unifying scaling diagram in the g' - q' 
plane as illustrated in Fig. [TU1 The flow lines clearly indi- 
cate that the phenomenon of "macroscopic charge quan- 
tization" is a universal feature of single electron devices 
that always appears in the limit where T goes to zero. 
Fig. [TU] furthermore displays all the super universal fea- 
tures of the 9 angle concept that previously arose in the 
theory of the quantum Hall effect. 

We end the paper with a conclusion in Section IVlII 



II. AES MODEL 



Action 



It is well understood by now that the AES model of 
the Coulomb blockade is a limiting case of the so-called 
universal theory of zero dimensional electron systems 
To start we briefly review the microscopic origins of this 
model. The experimental design of the SET is illustrated 
in Fig. The hamiltonian is split in three distinctly 
different parts 



Ti = Tig + Tl c 



E 

s—l.r 



H 



The first part is the free electron piece 



Tto = 



E 4 s) 4 s)t 4 s) 

k.s—l.r 



^ ^ d^ a da . 



(2) 



(3) 



The index s runs over of the reservoirs on the left hand 
side (I) and right hand side (r) of the island respectively. 
The subscript k denotes the electronic states in the reser- 
voirs and a those on the island. The ejr^, e a are the en- 
ergies relative to the Fermi level. 

The second term in Eq. (|2|) is the result of the Coulomb 
interaction between the electrons on the island 



H c = E c 




(4) 



E c = e 2 /(2(C; + C r + C g )) stands for the charging en- 
ergy and q = CgVg/e represents the external charge on 
theisland (see Fig. [Ud). 
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FIG. 1: a) Sketch of the SET device, b) Equivalent circuit of 
the SET. 
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The last part of Eq. describes the tunneling of 
electrons between the reservoir and the island 



H 



(s) 

The matrix t k ^ contains the amplitudes for tunneling 
between the reservoirs and the island. To characterize 
this tunneling it is convenient to introduce the following 
hermitcan matrices 



» 



w„wt 



dn 



h.c. 



(5) 



»(») 
9kk> 



4tt 2 



1/2 



S(ei s) )5(e^) (6) 



A: a:' 



(7) 



The first matrix acts in the Hilbert space of states of a 
single reservoir and the second one in the Hilbert space 
of states of the island. One should think of the delta- 
functions in Eqs. ©-0 as being smoothed out over a 
scale SE such that max{<5, 6^} < SE < T. Here, 8 
and (*>('' r ) stand for mean level spacing of single-particle 
states on the island and reservoirs respectively. 

The classical dimensionless conductance (in units 
e 2 /h) of the junction between a reservoir and the island 
can be expressed as follows^ 



/ j 9kk 

k 



(8) 



Therefore, each non-zero eigenvalue of g^ or g( 8 ) corre- 
sponds to the transmittance of some 'transport' channel 
between a reservoir and the island^ The effective num- 
ber of these 'transport' channels (N^) is given by 



1^ 9kk 

k 
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9 aa.' 9 a' a 



(9) 



The effective dimensionless conductance °f a 'trans- 
port' channel can be written as follows 



V a {s) a {s) V a [s> a (s > 

l^9kk'9k'k 9aa'9 a 'a 
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The dimensionless conductance g s then becomes 

a - a {s) N (s) 
9s — g c h Jv ch ■ 

In what follows we will always assume 

*r « i. 



(10) 



(11) 



(12) 



Notice that under these circumstances the conductances 
gi t r can still be large provided the effective number of 

channels N^b r ' ^ 1 is sufficiently large. 

We furthermore assume that the mean level spacing is 
negligible max{l,g},— and the charging energy 

is sufficiently large E c ^> 5 such that the effects of the 
exchange interaction can be ignoredi^i 



1. Path integral representation 

Given this sequence of limitations one can express the 
dynamics of the SET in terms of a single abelian phase 
$(t) with r standing for the imaginary time.-i This field 
generally describes the potential fluctuations on the is- 
land according to V(t) = i$(r). The quantum mechan- 
ical partition function Z can be written as a sum over 
winding numbers W according to 



Z[q]= e 2nqW Z 



(13) 



where Zyy is the integral over all paths 3>(r) that start 
with $(0) at r = and end with $(/?) = $(0) + 2ttW at 
r — (3 with (3 the inverse temperature: 



AY 



*(,3)=*(0)+27rW 



V$( T )e- Sd W- s *W. (14) 



Here, the action Sd describes the tunneling between the 
island and the reservoirs 



S d [§\ = f / d Tl dT 2 a(r 12 ) 
4 Jo 



,i$(Ti) — i*(T2) 



(15) 



where g = gi + g r and T12 = t\ — t 2 . The kernel a(r) = 
a(r + /3) in time and frequency representation is given by 

a(r) = -T 2 coscc 2 (7rTT) = - V |w rl |e-^" T (16) 



with w„ = 2nTn. The second term 5 C corresponds to the 
charging energy due to the Coulomb interaction between 
the electrons on the island 



(17) 



Finally, the exponential factor containing the external 
charge q in Eq. (p~3|) describes the coupling between the 
island and the gate of the SET. 



2. Functional integral representation 

An elegant formulation of the AES theory is obtained 
using the 0(2) field variable 



C(r) 



/ cos $ sin $ 
I sin $ — cos <£> 



(18) 



with Q 2 (t) = 1. The partition function can now be 
expressed as follows 



Z[q] = f V[Q]e 

JdV 



S[Q] 



(19) 



where the subscript dV indicates that the functional in- 
tegral is performed with periodic boundary conditions 



Q(0) = Q{P). 



(20) 
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The action is given by 



Instantons 



f 

S[Q] = / dndnifa) trQ(n) Q(r 2 ) 

Jo 



dr tr <7 V Qd T Q 



(21) 



with er denoting the Pauli matrices. The kernel 7(712) in 
frequency representation can be written as follows 



An 



1 

8E~ C 



(22) 



Alternatively one may express the action in terms of the 
0(2) vector field, <2(r) = a z (ti x (t) + icr y n y (T)), with 
n 2 (r) = 1. The integer valued topological charge of the 
system can be expressed in three different ways 



C[Q] 



4tt Jo 







dTtra y Qd T Q 



1 f 

J dre^rip ■ d T n v 



2tt Jo 



dr $. 



(23) 



This quantity is nothing but the number of times (W) 
the 0(2) vector field n is winding around. It is important 
to emphasize that both the periodicity statement of Eq. 
(|20| and the quantization of topological charge in Eq. 
(f23|) arc fundamental features of the AES theory that 
inherently describe the SET in thermal equilibrium. This 
theory only depends on the external charge q modulo 
2tt. If, for example, one splits q into a fractional piece 
— 7r < 6(q) < 7r and an integral piece k(q) (see Fig. [5]) 



%) 
2tt 



k(q) 



(24) 



then the quantum mechanical partition function only de- 
pends on the fractional piece 9(q) 



Z[q] = Z[6(q)/2n] 



(25) 



To extract the integral piece k(q) from the AES theory 
one must in general consider perturbing fields that take 
the SET out of thermal equilibrium. 



One of the most impressive features of the tunneling 
term of Eq. (p~5]) is that it possesses stable classical min- 
ima &w(t) for each topological sector W. We term these 
classical solutions "instantons" since they are completely 
analogous to Yang-Mills instantons^ The general ex- 
pression for <J>iy(r) is given b y 20 ' 21 



with 



ZlT 



z(r)z a 



-2-KiTT 



(26) 



(27) 



For instantons (W > 0) the complex parameters z a are 
all inside the unit circle and for anti-instantons (W < 0) 
they are outside. The classical action 



ir J 



fi*i 



(28) 



is finite and independent of the complex parameters z a 
which are the 2\W\ zero modes in the problem. 

In the limit where g — > 00 one may generally think 
in terms of a dilute gas of single instantons and anti- 
instantons. One identifies To = argzi/27rT as the po- 
sition (in time) of the single instanton whereas A = 
(1 — I 1 2 )/3 is the scale size or the duration of the po- 
tential pulse i$±i(t). The thermodynamic potential 
^inst = —T\nZ of the dilute instanton gas^ 2 - can be ex- 
pressed in a standard manner as an integral over To and 
A according to 12 ' 23 




12 3 4 

FIG. 2: Integer k(q) and fractional 6(q) pieces of q. 
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pn inst = - I dT I ^g(\)De-^ W+ ^ {T -^cos2irq. (29) 



r 



with D — 2er 1 ~ 1 and 7 w 0.577 the Euler constant. bation theory 2 ^ 
Here, the quantities 3(A) and E C (A) have the same radia- 
tive corrections as those obtained from ordinary pertur- _ ^ _ 2l n AA E (A) = E ( 1 — — In AA ) (30) 
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with A = gE c /w 2 D standing for the frequency or en- 
ergy scale. From Eq. (|30[) we obtain the renormalization 
group equations which to order g" 1 are given by 



P 3 = 



dg(\) 

d\n\ 



= -2- 



5(A) 



/3 C = 



dlnE c (X) 



din A g(X) 
(31) 

Here, we have included the perturbative contributio n 2 ^ 26 
of order l/g(X) into f3 g . Based on perturbation theory 
alone one expects that the quantum system is a good 
"conductor" at high temperatures 



5 (T) = -21n/?/£»l 
and an "insulator" at low temperatures 

ff (T)=cxp[-(/3/e) 2 l «1 



(32) 



(33) 



Here, z is a dynamical exponent that is as of yet un- 
known and £ denotes the dynamically generated correla- 
tion length in the time domain 



e = A-y 



(34) 



These standard ideas and expectations do not reveal 
much about the 9 angle concept on the strong coupling 
side, however. For example, there are the conflicting 
claims made by the semiclassical picture of the Coulomb 
blockade which say that the system displays a vanish- 
ing energy gap or a "quantum phase transition" when 
q passes through half integral values, see Section IIII Bl 
These conflicting scenarios raise fundamental questions 
about the exact meaning of the topological excitations in 
the problem and, in particular, the dilute gas of instan- 
tons written in Eq. (|29|) . 



III. PHYSICAL OBSERVABLES 

As pointed out many times in our previous work, tra- 
ditional instanton results such as Eq. (|29|) are of lim- 
ited significance since they merely describe the regular or 
non-critical pieces of the theory which are of secondary 
interest. In order to be able to understand the low energy 
dynamics of the SET and, in particular, the phenomenon 
of charge quantization one must develop an entirely dif- 
ferent approach to the AES model and reconsider the 
traditional renormalization group ideas in quantum field 
theory all together. 

Recall that conventionally one defines a renormalized 
theory by specifying how the ultraviolet singularity struc- 
ture of the bare theory can be absorbed in counter terms. 
There are many ways of doing this and normally, in the 
theory of critical exponent values in e expansions for ex- 
ample, one chooses a specific scheme based on computa- 
tional advantages. 

The infrared problems associated with the instanton 
angle 9 dramatically alter the physical objectives of the 
renormalization group. The extensive list of studies on 



the AES model is in many ways a reflection of what 
started many years ago in quantum field theory. There 
are the perturbative weak coupling analyzes, the instan- 
ton investigations as well as the various different at- 
tempts toward the strong coupling phase of the SET. 
Each of these distinctly different approaches to the AES 
model provide different pieces of knowledge in physics. 
They arc completely disconnected, however, and have 
physically very little in common. 

The basic idea pursued in the theory of the quantum 
Hall effect is to provide a unifying renormalization the- 
ory of the instanton angle 9 based on the response of 
the system to infinitesimal changes in the boundary con- 
ditions. This idea is very close to the criterion of An- 
derson localization originally proposed by Thouless^i It 
is also very close to 't Hooft's idea on duality based on 
twisted boundary conditions 2 ^ which states that gaplcss 
excitations must in general exist when 9 passes through 
odd multiples of n. Unlike these well known principles 
in physics, however, one now relates the sensitivity to 
boundary conditions to a set of "physical observables" 
that provide a very general definition of the renormal- 
ization behavior of the system. In the context of the 
quantum Hall effect these physical observables have pre- 
viously been recognized as the macroscopic conductance 
parameters of the system. 

The AES model is an interesting and highly non-trivial 
example where the theory of physical observables can be 
explored and investigated in great detail. In this Section 
we show that the problem of charge quantization in the 
SET is completely analogous to the robust quantization 
of the Hall conductance observed in the disordered elec- 
tron gas in two dimensions. It turns out that the AES 
model is extremely interesting in and of itself because of 
the the long ranged nature of the tunneling term or the 
non-local properties of the kernel 7(712) in Eq. (f2"Tj) . For 
the sake of simplicity we assume throughout the present 
Section that 7(712) is local in time and postpone the re- 
finements and extensions of the argument to Section [TV] 

As a trivial but very instructive example of our general 
definition of physical observables we study the isolated 
mesoscopic island in Section al! Bl This naive strong cou- 
pling example reveals much the conceptual structure of 
the instanton angle 9 and sets the stage for the remainder 
of this investigation. 



A. Background fields 



Consider a fixed background matrix field Uq(t) or 
Qo(t) — U a ' 1 a z Uo that varies slowly in time. We as- 
sume that <2o satisfies the classical equations of motion 
and carries a small fractional topological charge, i.e. Qo 
violates the boundary conditions of Eq. (|2TJ|) . The theory 
in the presence of the background field 



Z[q; Qo] = / V[Q]e 

JdV 



-S[U QUa 



(35) 
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then provides all the important information on the quan- 
tum system at low energies. To relate the background 
field action the appearance of an energy gap in the SET 
one must separate the constant pieces in Qq from the 
parts that couple to the the matrix field variable Q. Em- 
ploying the split of Eq. (j24|) and keeping in mind that 
C[Q] is quantized then one can write 

cxp {2mqC[U QUa 1 ]} = 

= cxp{27rik(q)C[Qo}+ie(q)C[U QUo 1 }} .(36) 

Using this identity one can split the theory of Eq. (155)) 
into pieces that are periodic and non-periodic in the ex- 
ternal charge q according to 

Z[q;Q ] = e 2 " fc W c [2o] Z[6(q)/2*;Qo\. (37) 

It is clear that only the periodic piece probes the sensi- 
tivity of the SET to changes in the boundary conditions. 
Provided the 7(712) is local in time one obtaines the ef- 
fective action in Qq in terms of a derivative expansion. 
The result is of the same form as the AES action itself 



SefflQo] 

S e [Q ] 



= In 



2mk(q)C[Q ] + S B [Qo] 
Z [%)/2tt; flo] 



Z [9(q)/27r] 



(38) 
(39) 



= J dr\ dr 2 j'(n2 ) tr Q (n) Q (r 2 ) 



- ie'C[Qn]+0{Ql). 

except that the bare quantities 7 and 9(q) are replaced 
by the effective expressions 7' and 9' respectively. As a 
criterion for a mass gap or energy gap in the SET one 
can now state that <Se[Qo] must vanish order by order in 
an expansion in powers of the derivative acting on Qq. 
This means that not only the 7' and 9' are exponentially 
small in (3 but also the infinite series of higher order terms 
not written in Eq. ([591 . Under these circumstances the 
effective action is given by 



Scff[Qo] = 2mk{q)C[Q \. 



(40) 



In the context of the disordered electron gas one identifies 
this result as the action of "massless chiral edge excita- 
tions." The quantity k(q) is recognized as the robustly 
quantized Hall conductance with sharp steps occurring 
at the center of the Landau bands (i.e. q = m + 1/2 with 
integer in). 

Presently, the background field Qq merely stands for a 
perturbing field that takes the SET out of thermal equi- 
librium. The quantity k(q), however, is identified as the 
robustly quantized quasi particle charge of SET. This 
quantity, as we shall see in Section IIV1 is in general very 
different from the averaged charge Q on the island. 



B. Isolated island 

To see these general statements at work we go back 
to the path integral representation of Section III A II and 



consider the simple problem with the tunneling conduc- 
tance g equal to zero. The classical equation of motion of 
the Coulomb term of Eq. ([T7|) is given by d 2 ^/dr 2 = 
which is simply solved by writing 



$(t) = 2irT(W - 



(41) 



The integer W generally stands for the integral topologi- 
cal sectors of the system and — l/2<0<l/2 denotes the 
perturbing background field with a fractional topological 
charge. We can write 



w ^ 

= Z[q]e- S ^W 



(3E C 



(42) 



The "effective action" 5 e ff[0] m Eq. (|42|) has the same 
general form as the original AES theory (in the absence 
of tunneling) 



«S eff [0] - -2mq'cj> - + 0{4?) 



(43) 



except that the bare parameters q and E c are now re- 
placed by the effective or "observable" ones q' and E' c 
respectively. It is readily seen that 



9 
1 

E 7 . 



1 din Z[q] 
C1 + 20E C dq ' 

_L (lA 1 d 2 In Z[q] 
E c { 2(iE c dq 2 



(44) 
(45) 



Similar expressions can be written down for the coeffi- 
cients of the higher order terms in iS e fi which in general 
arc irrelevant. 



1. Further evaluation 

To investigate the criterion for charge quantization 
written in Eq. (|4"0|) we must evaluate the observable the- 
ory of Eqs $2$ - (g5j) in the limit T = 0. Making use of 
the Poisson summation formula 



w 



2tt ^ v 



n) 



(46) 



one can express the partition function of Eq. (1421) as a 
rapidly converging sum over quantum numbers n accord- 
ing to 

Z\q-A\ = ^exp^Tm^-^n-g) 2 }. (47) 



We immediately recognize the grand partition function 
for Eq. (J4| with the integer n now standing for the num- 
ber of electrons on the island. The effective action can 
now be written as follows 

S off [cf\ = 2m{n)(j ) - 2ir 2 ((n 2 ) - (n) 2 )0 2 + O(0 3 ). (48) 
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FIG. 3: The thermodynamic potential of the isolated island 
at finite (dashed curve) and zero (solid curve) temperatures. 



Comparison with Eq. ([70)) shows that q' is none other 
than the averaged charge (n) on the island and E' c is 
related to the variance 



J = (n) 



1 



8((n 2 ) - (n) 2 ). 



(49) 



To obtain explicit expressions for q' and E' c we follow up 
on Eq. (J37J and split Eq. (|4"7[) into periodic and non- 
periodic parts in q according to 

Z[q- 0] = e 2 ***^ Z[d(q)/2n- 0] 



2tt 



^ exp <j 2vrmV - f3E c ( n' - 



(50) 



It is immediately clear that Z[8(q)/2w; <j)\ in the limit 
(3 — > oo is independent of </>. In complete accordance 
with the general statement of Eq. (|4TJ)) we conclude that 
the island for all values of — 7r < 6(q) < it develops an 
energy gap. The quantity q' or the averaged charge (n) 
on the island is quantized 



(n) = k(q) 



(51) 



with sharp steps occurring at half-integral values of q 
where the energy gap vanishes. The thermodynamic po- 
tential 



mi- 



(52) 



displays a "cusp" at half-integral values of q indicating 
that the transition is a first order one (see Fig. [3]). 



2. Renormalization 

This takes us to the most important part of this exer- 
cise which is to show that the physical observables gener- 
ally define the renormalization behavior of the island at 
finite T. Notice that Eq. ({50")) is dominated by the terms 
with n' = 0, ±1. Write 



Z[0(q)/2Tr; 0] = Z[e(q)/2n]e- 



4> 2 + O^ 3 ) 



(53) 
(54) 



then the explicit results for the thermodynamic poten- 
tial f2[g] and the physical observables 8' and E' c can be 
written as follows 



/3fl[q] = -lnZ[6(q)/2ir 
1(1- A 
7T dn[ q ] 



i(/3£ c )(l-A ) 2 -ln 



1 + e -(^) A « 



0(g)- 
± 



E c dq 
2n 



e (pE c )A 
1 



(3E' C (3E C 



2tt 



1 



1 - 



2E C 
2^ 



dq 2 J 



(55) 



(56) 



(57) 



Here, ± denotes the sign of 9{q) and A is recognized 
as the the dimensionless energy gap of the island which 
vanishes near the critical point according to 



A 



1 - 



(58) 



Finally, we express the response quantity 6' in differential 
form and obtain (see Fig. [4]) 



de 1 e 1 

rfln/3 _ 2^ 



[2tt - |6>'|] In 



2tt- 



(59) 



This result clearly translates the physics of the isolated 
island in the language of the renormalization group. No- 
tice that the quantity E' c in Eq. (|57[) does not lead to 
more complex renormalization behavior since it is ex- 
pressed in terms of 6' alone. The same is true for the 
higher terms in Eq. ([5^| . 

We identify two different kinds of strong coupling fixed 
points, a stable one at 8' = and a critical one at 6' = 

±7T. 

1. Near the critical fixed point 6' = ±7r we find 

p e (9') =±tt + 6' (60) 

which is a standard result for a first order transition 
in one dimension. Eq. (|60p determines the energy 
gap exponent v according to 



1 
p 



d(3 e 



89' 



= 1. 



(61) 



Eqs (|4"B"]) and ([57]) tell us that near criticality the 
charge on the island is broadly distributed, i.e. the 
fluctuations are of the same order of magnitude as 
the averaged value 8'. 

2. Near the stable fixed point at 8' = we find 

Pe(0')=8']n\8'\ (62) 
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FIG. 4: Pe function (Eq. (|59j) ) for the isolated island. 



indicating that the averaged charge q' on the iso- 
lated island is robustly quantized with corrections 
that are exponentially small in (3, i.e. 



q' = k(q) + — = k{q) ± e"^ A ° . (63) 

Z7T 



Similarly, the root-mean-square fluctuations in 9' 
as well as the higher order moments all render ex- 
ponentially small in (3. 



IV. GENERAL RESPONSE THEORY 

Armed with the insights obtained from the isolated 
island we next address the AES theory with finite values 
of g. To discuss the tunneling term iS ( j[ ( I > ] with varying 
boundary conditions on the $ field one must generalize 
the expression for the kernel a(r 12 ) in Eq. p^|) which is 
periodic in time. Write 



o*(n 8 ) = -Y\ e~ i ^ +2 ^ n ' K + 2tt1>| (64) 

7T * — ' 



with —1/2 < <fi < 1/2. The appropriate result for the 
tunneling term S e i in Eq. (|16p is then obtained if one 
replaces a(ri2) by the following expression 

a(r 12 ) -> e^ T ^a^(r 12 )e-^ T ^ (65) 
= a(T 12 ) + 2T 2 |0|-2iT 2 0cot(7rTT 12 ). 

Eq. (f6"5")) essentially tells us that one cannot insert a back- 
ground field Sd[®] — > Sd[$ + $o] with $ = 2vrT(/i carry- 
ing a fractional topological charge unless one changes the 
kernel a(ri 2 ) into ol${t\2)- Given Eq. ()65[) it is straight- 
forward to discuss the effect of the more general back- 
ground field 



and the result can be written as follows 

P 



S d [$ + $ ] = f j dT!dT 2 e i *(^)-i*(T2)+i*o(r 1 )-i4 (r 2 ) 


P 

xa (r 12 ) = £-[ dndrze^-^ 



o 



(67) 



Notice that Eq. (j6"B"|) now satisfies the classical equations 
of motion of the AES theory as a whole, i.e. not for only 
the isolated island as discussed in the previous Section 
but also for the theory in the presence of tunneling. We 
will next embark on the distinctly different ways of han- 
dling the background field methodology depending on the 
topological charge of the field <f>o- 



3. <J>o with fractional topological charge 

By taking <j> = 2irT<frT or ui m = then Eq. 
can directly be used to probe the sensitivity of the SET 
to changes in the boundary conditions. Introducing the 
two-point correlation function 



rP rP 

D(iw„) =T / dridr 2 e^" T12 
Jo Jo 



-i$(Ti)+i$(T 2 ) 



(68) 



then to lowest orders in the 4> we obtain the total effective 
action 



(69) 



with <S c ff[<^] given by Eq. (|48[) and below. Keeping in 
mind that — !/2<0<I/2we split the sum in Eq. (|S"9"|) 
in n = and n ^ parts and we immediately obtain 



«5 to t [4>] 



\<j)\-2mq'<j) + 5S to t[4>] 



(70) 



where <5<5>tot stands for all the higher order terms in (j. 
The physical observables g' and q' are given as follows 



g' = gTD(iO), 

q = Q - —Ty^ ]mD(iu) n ). 

Z7T ' — ' 



ri>0 



Here we have introduced the quantity 



Q = q 



2E C 



= q 



1 dn[ q ] 

2E C dq 



(71) 
(72) 



(73) 



*o(r) = (Lu m + 2irT(/))T 



(66) 



which generally stands for the averaged charge on the 
island. We see that in the presence of tunneling the av- 
eraged charge Q is different from q' which we now identify 
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with the quasi particle charge of the SET. We expect that 
the new quantity q' is quantized and, along with that, the 
quantity g' as well as all the higher dimensional terms in 
(5<Stot render exponentially small in the limit where (3 goes 
to infinity. 



4- Higher dimensional terms 

To obtain the leading order corrections in <5<Stot one 
needs the four point correlation function 



D(iuj n ,iuj m ) = T~ 



12 J 34 



(74) 



The effective action SStot up to the third order in cf> can 
be written as follows 

^totM = ^ 2 + + 0(d) 3 ) (75) 



where 



1 



1 9 / , 



0E' C (3E C 8q 



(V-Q) +^-D(i0,i0) (76) 



g 2 T 2 



8ir 2 

1 1 Og' ig 2 T 2 

$F' C ~ 4(3E C ~dq + 8tt 



^ sgn{u> n uj m )D(iu! n ,iuj m ), 



"Y] sgn u> n D(iu n , i0).(77) 



n#0 



A more detailed discussion of the higher order terms 
will be presented elsewhere^ 



5. $o wif/i integral topological charge 

A less obvious way of probing the energy gap in the 
SET is obtained by putting <j) = and, instead, we 
consider background fields with an integral topological 
charge only, i.e. <&o = u m r. Notice that this choice of 
$o is a special case of the instanton solution of Eq. (|26|) 
with W = m but with all the parameters z a put equal to 
zero. Even though this background field $ can formally 
be absorbed in a redefinition of the $ field one can nev- 
ertheless proceed and define the effective action <5>tot[$o] 
by expanding in powers of $o or uj m . Provided one finds 
a way to analytically continue the discrete Matsubara 
frequencies to fractional or infinitesimal values the final 
results are again a measure for the sensitivity of the SET 
to changes in the boundary conditions. 

To start we consider the effective action at a tree level 



m 2 . (78) 



S[$o] = ^ M ~ Zniqm + — 

We expect that the exact result retains the general form 
of Eq. (|75|) except that the bare parameters <?, q and E c 



are replaced by effective or observable ones. To lowest 
orders in m one can write 

StotM = S eS [m]+K(iw m )-K(iO) (79) 

where <S c ff [m] is the same as Eq. (|48p with <fi replaced by 
m. We have introduced the quantity 



K(iw n ) 



9_ 

'4/3 



d,Tid,T2 e 1 ' 



3 l*(Tl) — i$(T2) 



(80) 

To expand this theory in terms of a series in powers of 
Lj m we make use of the analytic properties of response 
functions. Specifically, following the standard prescrip- 
tion iuj m — > lj + i0 + we analytically continue the discrete 
set of imaginary frequencies iuj m in Eq. (|80[) to real ones 
uj and subsequently we can take the limit u> — > 0, see 
Section TlV Al The following total result is obtained for 
the effective action up to order m 3 



\m\ — 2iriq'm 



PE' C 



2ni 



(81) 



Here, the quantities q' and g' arc given in terms of Kubo- 
like expressions as follows 



An Im 



Q + Re 



dK R {u) 
dK R {uo) 



du> 



(82) 
(83) 



K r (lu) denotes the analytic continuation of K(iuj m ) and 
the limit uj — > is understood. As before the Q de- 
notes the averaged charge on the island (see Eq. (J73J)). 
Similar results can be written down for E' c as well as F' c 
but the expressions are considerably more complex than 
those of Eqs (|76[) - ([77]) and the details will be presented 
elsewhere. 

Even though the physical observables of this Section 
are formally different from those in the preceding Section 
they should nevertheless define the same renormalization 
behavior of the SET. In particular, in the presence of an 
energy gap the physical observables in both Eq. ([70)) 
and Eq. (|8Tj) should all scale to zero as /? goes to infin- 
ity except for the quantity q' that can take on arbitrary 
integral values. The main advantage of Eqs ([70]) - (|83[) . 
however, is that they directly lay the bridge between the 
background field methodology on the one hand, and re- 
sults obtained from ordinary linear response theory on 
the other, see Section TlVBI 



A. More about response functions 

The function K r (lu) can elegantly be expressed in 
terms of the retarded propagator D r {lu) which is the an- 
alytic continuation of D(iu> n ) in Eq. (|68p . In Appendix |A1 
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we derive the following relation 



and the corresponding Matsubara correlation function by 



j dexde 2 _ n b (e 2 ) - n b (ei) lmD R^ 



K"{u>) = g 



e 2 

47T 3 ei - e 2 + w + i0+ 



(84) 



where rc.(,(e) = [exp(/3e) — l] -1 denotes the Bose-Einstein 
distribution. A detailed computation of K R (u>) in the 
weak and strong coupling regimes is presented in Ap- 
pendix [D] . 

Given the function K R (ui) one can obtain the response 
parameters g' and q' from Eqs. (j82")l and {53]). However, 
it is more convenient to express these quantities directly 
in terms of D r {uj) according to 



de dn b R 
— e— lmD H e , 

7T 0€ 



de den b 
4^2 ~dT 



ReD R (e). 



(85) 



(86) 



B. Linear response 

In this Section we establish the contact between the 
response quantities g' and q' and the well known ex- 
pressions for the SET conductance G and the non- 
symmetrized current noise Si. 



K s (iiu n ) = / dTe^ T (T T X<- s \T)X^(0)). (91) 
Jo 

Repeating the same steps that led to the AES action 
starting from the Hamiltonian of Eq. ((2)) we obtain 



K s (iu) n ) 



4/3 



f 

/ d Tl dT 2 e w ^a{T 12 )D{T 21 ). (92) 
Jo 



Comparison with Eq. (|80|) yields K s (iui n ) 
(g s /g)K(iu n ) or, equivalently, 

K R {u) = (g B /g)K R (w). 



(93) 



Based on the continuity equations for the current / = 
Ii = —I r = GV we finally find the SET conductance G 
in units of [e 2 /h] according to 16 ' 31 i 32 



G 



9i 9r , 



(gi + gr) 2 ' 



(94) 



with g 1 given by Eq. (f8"2"|) or ([55]) . Therefore, except for 
the constant gig r /(gi + ffr) 2 the conductance G is none 
other than the observable g' that measures the sensitivity 
of the SET to changes in the boundary conditions. 



2. Quantum current noise 



Similar to Eq. (|89|) we obtain the real part of the re- 
tarded correlation function as follows 



1. The SET conductance 

If one applies a voltage difference V = V r — Vi between 
the reservoirs then the tunneling part TL^^ of Eq. ^ 
becomes time dependent^ 



W =1 W e -ieKt +Jf Wt e ieKi ) 



(87) 



The operator for the current I s that flows from a reservoir 
to the island can be expressed as follows 



l s = ej t Y, 4 s)t «i s) = -ieX^e-^ + h.c. 



To the lowest order in 1/Nj$ we find 

la = -I [ dt'([lW(t),H^(t')}) 
J — oo 

= -2elmK R (-eV s ). (89) 



The retarded correlation function is given by 

/•OO 

K R (u) = i dte iut ([X^(t),X^(0)]) (90) 
Jo 



ReK R (-eV a ) = — dt'{[X s (t),X s (t')}). (95) 



2e 2 



The quantity q' in Eq. (fBU)) or (|86[) can therefore be ex- 
pressed in terms of the current-current correlation func- 
tion^ 

ql = Q _ ^m + grf d r° dt{[ma{t)]) (96) 

2#Zffr OV J_ oc 

in the limit where V goes to zero. We have thus 
found a novel interpretation of the so-called antisym- 
metric current-current correlation function that in dif- 
ferent physical contexts has attracted a considerable 
amount of interest over the years^l. Introducing the non- 
symmetrized current noise^ 



Si(lo,V) 



then one can also write 



dte- Mt {l(t)I(0)) (97) 



gig r 



2tt uj dV 



Here, PV denotes the principal value and the limit V — > 
is understood. Eqs (|M|) and (f9"B")) arc amongst the most 
significant results of this investigation. 
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V. WEAK COUPLING REGIME, g' > 1 
A. Perturbation theory 

At a gaussian level the AES action in frequency repre- 
sentation is given by 



n>0 



2ir 2 T 



n 2 $„$_ n . 



(99) 



To lowest order in an expansion in g the following result 
for D(iuj n ) is obtained 



D(iu n ) = 13 



g £g a + 2-K 2 Ts 2 /{gE c ) 



1 



1 



2iri 

g \i\uj n \ i\Lu n \+igE c /i: 



(100) 



Eq. (|29[) . This leads to the much simpler expressio: 



/?fi inst = -4/^ce- fl / 2 ln ; 



2ir 2 et 



cos 2irq. 



(104) 



With the help of Eq. ([73)) we immediately find the tem- 
perature dependence of the average charge on the island 
and the result is 



Q(T) =q - — e~ s/2 ln 



2-n 2 e<T 



sin27rg. 



(105) 



To find the quantities q' and g', however, we still have 
to evaluate the instanton contribution to the correlation 
function D(iuj n ). For this purpose we first consider the 
expectation of an arbitrary operator O which can be ex- 
panded to lowest order in the topological sectors W = ±1 
according to 



Using the representation 5 n o = hm r)(i(j n +17) 1 one can 

perform analytic continuation to real frequencies and the 
retarded correlation function becomes 



(O) 



1 



Z[q] 



(106) 



D R (u) = P 



2ni 



g 2tt 2 T 
1 



n 



lim - 

»?^o ijj + 77 + i0 + 

1 



£J + i0 + lu + igE c /n 



(101) 



Having carried out integration in Eqs. ([55]) and (fB"6"|) with 
the help of the identity: 



where 



O w = I D$(t) 0($) e 

/$(^)=$(0)+27rW 



-S«t[*]-5 C [*] 



(107) 

Similarly, we expand the partition function according to 



o x 2 + n 2 z 2 sinh 2 x 2\z 



= —-1 + 1^(1. 



(102) 



where ip{z) denotes the Euler di-gamma function, we ob- 
tain 



nE p 7+1 

</(T)= 5 -2m^£_ q'(T) 



(103) 



Here, 7 = —"0(1) ~ 0.577 denotes the Euler constant. 
The result for g' was originally obtained in Ref. [24| more 
than two decades ago. The quantity q' , on the other 
hand, is unaffected by the quantum fluctuations to any 
order in an expansion in powers of 1/g. To establish the 
renormalization of q' (0-renormalization) it is necessary 
to include the non-perturbative effects of instantons. 



Z[q]^Z (l + e 2 ^^ + e- 2 ^^ 



(108) 



Eq. (| 1 06[> can therefore be split in a W — part and an 
instanton part 



(O) « (O) + (0) inst . 



(109) 



Here, {O) = O /Z and 



^0 



(110) 



B. Instantons 

Since the infrared of the dilute instanton gas is well 
defined one can proceed and evaluate the integrals in 



In the semi classical evaluation of Eq. (|110p it suffices to 
replace the operator 0[$] in the integrand of Eq. <\W7\i 
by its classical value 0[$vH- Tnc result for Eq. (|110p 
can then be written in the typical instanton form 
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rP rP j\ 

(0) ins t = V / dr / — - (O) ] fl(A)D exp 



4 9W + «A) 



(111) 



r 



where ©[$±1] generally depends on the position To and 
scale size A of the instanton/anti-instanton. 

We next apply these general results to the correlation 
function D(iuj n ). The operator specific parts of Eq. 
are computed to be, to the leading order in 1/g, 



dr [0[$ w ] - (O) ] 



Sn,0 




with denoting the Heaviside step function. Inserting 
this result in Eq. (|29|) and performing the integral over 
A we find the following result for the instanton part 



AnstC^ 



9 2 E C _ 
7r 2 T 2r 



ff/2 



X (1 - $n,o) ~ 



5 nfi cos2Trq-TriTe 27Ttqssnn 
1 1 



i\uj n \ + 2niT 



(113) 



Performing analytic continuation to real frequencies we 
obtain 



AnstM 



9 2 E C 

7T 2 T 2 ' 



-9/2 



cos 2nq lim 



n 



niTe 



1 



1 



lu + i0+ uj + i2ttT 



(114) 



Using Eqs. and (|86[) we find the following non- 

perturbative corrections to g' and q' 



g 3 E c 



6T 



-9/2 



COS 27T(7, 



g 3 A 

24ttT 



9/2 



sin27rg. 



(115) 
(116) 



Here, the expression for Q(T) is given by Eq. (|105p . 

Combining the perturbativc and non-pcrturbative con- 
tributions of Eqs. ([TUB")) . pUS]) . PIS)) and (fTTrJ)) we obtain 
the final total result for the temperature dependence of 
g' and q' 



gE c e 



7+1 



9 3 E C 



q'(T)=q- 



g 3 E c 

24irT 



9/2 cos27r(7, (117) 
9/2 sin 2irq. 



(118) 



Several remarks are in order. First of all, we notice that 
the amplitude of the the oscillations in q' with varying q 
are much larger than those in the averaged charge Q(T). 
Eqs (|117[) and (|118[) arc generally valid in the weak cou- 
pling phase of the SET g < 1 such that T > g 3 E c e~ g / 2 . 
The results are completely analogous to the instanton 
corrections to the conductances a' xx and <j' xy in the the- 
ory of the quantum Hall effec t 37 ' 38 that have recently 
been investigated experimentally^ It should be men- 
tioned that Eq. ()115|) coincides with the earlier compu- 
tations reported in Ref. [Til ]. 



renormalization 



To leading order in 1/g one can express Eqs. (|117p and 
(|118[) in the following manner 



g'(T) = g{T) - Dg 2 (T) e - g{T)/2 cos2tt 9 
D 



q'(T) =q- ^g 2 (T)e- 9{T)/2 sm2irq. 

4-7T 



Here, D = (ir 2 ^e" 7 " 1 
and 

9(T) = 



(119) 
(120) 

0.68 is a numerical constant 



2 In 



gE c 

6DT 



(121) 



contains the perturbative quantum corrections to leading 
order in 1/g. It is important to emphasize that same re- 
sults of Eqs. (|119|) and (jl20[) arc obtained if one employs 
the much simpler expressions for g' and q' defined in Eqs 
(fTTj) and (f72")) . The only difference is the numerical value 
of D which now equals D = 2 exp (—7). At the same 
time, the charging energy E c in Eq. (|121[) is replaced by 
(6A 2 )£ C . 

Expressing Eqs. (|119[) and p20[) in differential form 



Pi 



4f~ =-2--- Dq' 2 e- 9 '/ 2 cos 2nq' (122) 



dq' 
din /3 



_ Rg'2 e -9'/2 sin27r(? '(l23) 

47T 



we obtain the renormalization group functions (3 g q — 
(3g,q(g',q') of the AES theory on the weak coupling side. 
We have included the two loop correction 2 ^ in the per- 
turbative part of Eq. ([1221) . 

Eqs (|122|) and ()123|) arc amongst the most important 
results of this investigation. The results clearly demon- 
strate that instantons are the fundamental topological 
objects of the AES theory that describe the cross-over 
behavior of the SET between the conducting phase at 
high temperatures and the Coulomb blockade phase that 
generally appears at much lower temperatures only. 
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VI. THE STRONG COUPLING PROBLEM, 

<?'<! 

A. Effective action for 8 ~ -n 

For small values of the tunneling conductance g we can 
simplify the Hamiltonian of Eq. ((2]) near the degeneracy 
point 9 = 7r or q = 1/2 by employing a projection onto 
the states with Q = k(q) and Q = k(q) + 1 of the isolated 
island^ The projected Hamiltonian can be written as 
follows 



Ti. = TCo + T~tc + T~C 



(0 



where 



E c (k-qf + - 



?4 s) = E4 s 24 s)t ^- 

An 



■4» 

+ h.c. 



(124) 

(125) 
(126) 




-o 



FIG. 5: The pseudofermion self-energy: solid line denotes 
G(ie n ) whereas wavy line stands for the interaction a(iu>„) 
(see text). 



is similar to the XY case of the Bose-Kondo model 
for spin s — l/2i 43 i 44 i 45 Notice that the spin operators 
4>(t)o-±i[>(t) in Eq. Q127P are the same as the AES oper- 
ators exp(±z<I>(T)) projected onto the states Q = k and 
Q = k + 1 of the isolated island. The projection onto the 
Hamiltonian of Eq. (|124[) is justified as long as g <C 1, 
\q - k - 1/2| < 1 and (3E C > lM 



Here, A = E c (l - 9/w) > 0, S denotes the spin s = 1/2 
operators and S ± — S x ± iS v . 

A convenient representation is obtained by using 
Abrikosov's two-component pscudo fcrmion fields ip and 
^i 40 i 41 After integration over the electronic degrees of 
freedom one arrives at the following effective action to 



leading order in l/N ( 



is) 
ch 



B. Leading logarithmic approximation 



In what follows we limit the analysis of Eq. Q127P to the 
so-called leading logarithmic approximation. This corre- 
sponds to the one-loop renormalization group procedure 

of Rcfs. HE!. 



^2" + / dT ^'[ d r-V+J<?z 



4 7o 



dn dr 2 a ( n 2 ) [ ( T\ ) a _ -0 ( n ) ] ( t 2 ) cr + V» ( r 2 ) ] . 



(127) 



Here, Oj with j = x, y, z stand for the Pauli matrices and 
a± = (a x ± ia y )/2. We have introduced the chemical 
potential rj such that the limit rj — > — oo is taken at the 
end of all calculations. This procedure ensures that only 
the physical states with pseudo fermion number N p f = 1 
contribute to the quantities of physical interest. Follow- 
ing the prescriptio n 40 ^ 1 



Z = lim 



d 



Z 



(128) 



we obtain the physical partition function Z from the 
pscudo fcrmion theory Z p f. Similarly, we extract the 
physical expectation (O) according to 



(O) 



lim 

Tj— > — OO 



d 



Z de^ 



(0) P f + (0) P f 



(129) 



The brackets (. . . ) p f denote the average with respect to 
the theory of Eq. fT2T]) . 

In what follows we employ the effective action of 
Eq. (|127p to investigate the phenomenon of macroscopic 
charge quantization as well as the renormalization be- 
havior of the SET on the strong coupling side. Eq. (|127p 



1. Pseudofermion Green function renormalization 



Using Eq. (|127|) we find the following expression for 
the pseudofermion Green function for g = 



G ±(ie n ) = ie n + r]T y, 



(130) 



where e n = irT(2n + 1). The pseudofermion Green func- 
tion can be expressed in terms of the self-energy Sj- 

G^}(ie n ) = ie n + i] =p y - S±(ie„). (131) 
It is convenient to parameterize the self-energy as follows 



S±(ie„) = (ie n + V )[l - 7 («e„)] T [1 - 7a(^n)]y (132) 
The pseudofermion Green function now becomes 

G± 1 (ie„) = (i^ + vh( ie n) T 7s(* e n)y- (133) 

The leading logarithmic approximation corresponds to 
the simplest diagram for the self-energy shown in Fig. [5] 
This leads to the following equation 



47r ^ iuji + ie n + f] ± — - S T (iw; + ie„) 

(134) 
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which has to be solved self consistently. Recall that there 
is no renormalization of the interaction line (see Fig. ([5])) 
because of the absence of closed fcrmion loops in the 
pseudofermion diagrammatic technique, i.e. their contri- 
bution vanishes in the limit r/ — > —ooi 40 i 41 

With logarithmic accuracy we see that both 
7 and 7s depend on the single variable x = 
In A/ max{A7 s /7, \ie n + rj\\ where A is an arbitrary high 
energy cut-off. Then from Eq. (|134[) we obtain 



7(z) = H 

Js(x) = 1 - 



9 

' 4tt 2 

.9 

4tt 2 







dy 
dyisjy) 



and, hence, 



1 ( x ) = 1 -\ x ) = (i + JL x ) 



1/2 



(135) 
(136) 

(137) 



2. The partition function and the average charge 



Using Eq. (fl"2"5|) we find 

z = e J3E c (k- q f e PA/2 Um e -fh, 
rj—> — oo 



e„,<r=± 



G a {ie n ). 
(138) 

Given the vertex functions j(x) and "f s (x) it is now trivial 
to evaluate the partition function 



2 



Here. 



7= (1+ 9 m A 

2tt 2 max{A',T} 



1/2 



(139) 



(140) 



and A' = A/7 2 stands for the renormalized energy gap 
between the ground state and first excited state. 

The average charge on the island is expressed in terms 
of the magnetization M = {S z ) of our spin model 



Q{T) = k+\ 



where 



M 



M{T) = _ tanh ^_, 



(141) 



(142) 



Eq. (|142[) has originally been obtained in Ref. [lg] using 
slightly different techniques. Evaluating the result at T = 
we find 



Q(T = 0) = k 



1 



1 



1 + 



(143) 



which is the familiar result of Matveev^ It does not 
resemble the simple expression for the averaged charge 
on an isolated island. This charge is, in fact, no longer 
quantized when the tunneling conductance g is finite, see 
Section rvTCTl 




FIG. 6: The two-point correlation function: solid line de- 
notes G(ie n ) whereas black triangle denotes the vertex func- 
tion T(ie n ) (see text). 




FIG. 7: The first non-trivial contribution to the vertex func- 
tion T(ie n ) (see text). 



3. The correlation function D r (lj) 

The diagram for the two point correlation function 
D(iuj n ) is shown in Fig. [51 Because of the peculiar form 
of the pseudofermion interaction which couples the <r_ 
and (7+ it is readily seen that the lowest order contri- 
bution to the vertex function T(ie n ) is proportional to 
g 2 ln(A/ max{T, A'}), see Fig.H Within the leading log- 
arithmic approximation one can therefore put T(ie n ) = 1 
and, hence, 



D{iui n ) 



xG + (ie m + iuj n ) 



2cosh/3A'/2 deP* 

tanh/3A'/2 1 



(144) 



7 2 iu> n - A' ' 
After analytic continuation to real frequencies we obtain 



tanh/3A72 1 
i H f,,A — — — L 

7 2 ui - A' + i0+ ' 



D R (u) = - 



(145) 



C. Physical observables 



Given Eq. (|145[) . it is possible to evaluate integrals in 
Eqs. CLE]), dZU) (definitions I) and in Eqs. (85]) and f86]) 
(definitions II) and compute the response functions in 
the strong coupling phase. For both definitions (I) and 
(II) we obtain the same expression for the quasi particle 
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charge q' 



I -i 2 3 A' 
q- = y + —J-tanh^ (146) 



k(q) 



2 7 2 2 
1 



oJA> + ! ■ 



(147) 



However, different expressions for the response quantity 
g' are obtained and the result is 



fJi 
9'n 



g tanh/3A'/2 



2 7 2 
9 



8A'/2 
3A' 



2 7 2 sinh 3 A 1 



(148) 
(149) 



Eq. ()149|) coincides with the result found in Ref. [l6( us- 
ing a somewhat different approach. It furthermore cor- 
responds to the sequential tunneling approximation of 
Ref. [11]. 

There are several interesting conclusions that one can 
draw from these findings. First of all, we see that as T 
approaches absolute zero Eqs (|146|) - (|149[) are indepen- 
dent of g and precisely coincide with the results obtained 
for the isolated island. Unlike Eq. (|143[) . for example, we 
find that the new quantity q' is robustly quantized with 
infinitely sharp steps occurring when the external charge 
q passes through k + 1/2. 

However, Eqs (| 148)) and (|149[) for the response quan- 
tity g' do not unequivocally predict an exponential de- 
pendence on T when 3 A' S> 1. Moreover, as shown in 
Appendix [C] the corrections to the quantities g' n and q\ 
to second order in g do not demonstrate an exponential 
dependence on T when T vanishes. 

This means that the strong coupling expansion in pow- 
ers of g generally does not provide access to the Coulomb 
blockade phase where the SET develops an energy gap. 
The validity of the leading logarithmic approximation is 
therefore limited to the quantum critical phase 3 A' < 1 
which for our purposes is the most significant regime of 
the SET. 

This takes us to the most important part of this ex- 
ercise which is to employ Eqs (|146l) - (|149l) in order to 
extract the scaling behavior of the SET on the strong cou- 
pling side. Since our physical observables are essentially 
defined for finite size systems (i.e. finite 3) they should 
in general be distinguished from the ordinary thermody- 
namic quantities of the quantum spin system that are 
normally being considered. Emerging from Eqs (|146D - 
(|149p there are two distinctly different renormalization 
group schemes, to be discussed further below, that pro- 
vide complementary information on the quantum system 
at zero temperature and finite temperatures respectively. 



RG at zero temperature 



Eqs (|146p - (|149p clearly show that two renormaliza- 
tions are in general necessary to absorb the ultraviolet or 



high energy singularity structure of our spin system, i.e. 
one renormalization associated with the coupling con- 
stant (tunneling conductance) g and one associated with 
the "magnetic field" (energy gap) A. From the expres- 
sions at zero T 



9_ 

V 



9 



1 



-Mn^- 

2tt 2 111 A' 



A' 



A 



Y 



1 + -K In ±- 



(150) 



we obtain the following renormalization group 3 and 7 
functions to one loop order 



0a 



dg 
din A 



9 

2tt 2 



din A 



7A 



din A 2tt 



(151) 



Employing the method of characteristics one can cast the 
thermodynamic quantities of the quantum spin system 
at T — in a general scaling form. For example, the 
magnetization M with varying "magnetic field" A can 
be expressed as follows 



M(A) = M o /(AM o 0. 



(152) 



Here, the functions M$ and £ with varying g are deter- 
mined by the 3 and 7 functions according to 



d_ 8 

'8A +Pg dg 
d_ 
' dg 



lA + d g — ) M = 0. 



(153) 
(154) 



One finds, for example, that the characteristic time scale 
£ of the SET is given by 



(155) 



which has the same meaning as the weak coupling expres- 
sion of Eq. ([34]) . The quantity Mo and the scaling func- 
tion f(X) within the one loop approximation are given 
by 



Af = l/2.g, /(A) = 27r 2 ln- 1 A- ] 



(156) 



The result essentially tells us that the spontaneous mag- 
netization only exists for the theory with g = but it 
vanishes for any finite value of g. In terms of the AES 
model this means that the averaged charge Q on the is- 
land is no longer quantized when finite values of the tun- 
neling conductance are taken into account. 



2. RG at finite temperatures 

We next specialize to the physical observables at finite 
temperature. Since the expressions of Eqs (|14€>|) - Q149P 
are universal for 3 A' -C 1 



^^(l + ^ln/M 



(157) 
(158) 
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FIG. 8: The function j3 g with varying g' along the critical line 
q — k + 1/2. An interpolation between the weak and strong 
coupling branches has been drawn as a guide for the eye, see 
text. 



FIG. 9: The function df3 q /dq' with varying g' along the crit- 
ical line q' — k + 1/2. An interpolation between the weak and 
strong coupling branches has been drawn as a guide for the 
eye, see text. 



we immediately obtain the finite temperature (3 functions 
along the critical lines q' = k + 1/2 for g' -C 1 according 
to 



P 3 = 
P q = 



dg' 
d\n(3 

dq' 
din (3 



9 



(159) 
(160) 



These results should be compared with Eq. (|59"|) obtained 
for an isolated island. We see that the critical fixed point 
of an isolated island is the critical fixed point of the AES 
theory as a whole with the SET conductance g' now play- 
ing the role of a marginally irrelevant scaling variable. 

Next we compare Eqs (|159[) and (|160|) with the the 
weak coupling results of Eqs (fl~22|) and (fl~23)) . In Figs [8] 
and[5]we plot the functions (3 g and d(3 q /dq along the crit- 
ical line q' = k + 1/2. A simple interpolation between the 
weak and strong coupling branches indicates that both 
these functions decrease monotonically as g' increases. 

Finally, it is not difficult to understand why the 
Coulomb blockade phase of the SET is beyond the scope 
of the present investigation. For example, given the fact 
that the theory with q' ss k and g' 1 develops an en- 
ergy gap then one generally expects (see also Eq. ([62]) ) 



A, = f/ ^ g' 

f3 q = (q' -k)\n\q' -k\ 



(161) 
(162) 



which cannot be obtained using ordinary perturbation 
theory in g' . 



VII. SUMMARY AND CONCLUSIONS 

To summarize the results of this investigation we have 
sketched, in Fig. [101 a unifying scaling diagram of the 
SET in the g 1 - q' plane. This diagram is based on the 
strong coupling results of Eqs ([59")) , (|159[) and (|160[) and 
the weak coupling results of Eqs (|122p and (|123p . 

The universal features of this diagram are the quan- 
tum critical fixed points located at q' = k + 1/2, g' = 0, 



and the stable fixed points at q' = k, g' = that describe 
the "macroscopic charge quantization" of the SET. The 
results are in accordance with the concept of super uni- 
versality that has previously been proposed in the context 
of the quantum Hall effect.— 

We have established the relation between the quantity 
g' and the ordinary SET conductance G that one nor- 
mally obtains from linear response theory. The quantity 
q' is new and can similarly be expressed in terms of the 
antisymmetric current-current correlation function. 

The quantization of q' is an interesting and important 
challenge for experimental research on single electron de- 
vices. There are, however, other ways of experimentally 
probing the quasi particle charge q' of the SET. In Sec- 
tion [VlTX] below we will summarize the quantum critical 
properties of q' and point out how they arc directly mea- 
surable in the experiment. 

We conclude this paper with Section lVII Bl bclow where 
we discuss in some detail the physical mechanism that is 
responsible for changing the the quasi particle charge q' 
of the SET as q passes through the critical point. 



A. Quantum criticality 



Eq. (|157p is the maximum value of g'(T) as one varies 
the value of q. This maximum value vanishes logarithmi- 
cally in T according to 



9max 

(T) 



9 
2 7 2 



_g_ 

2tt 



2 ln/3A 



= tt 2 !!!- 1 <c 1. (163) 



Similarly, Eq. (|158[) determines the maximum slope of 
the quasi particle charge q'(T) with varying q. This slope 
diverges according to 



7'(T) 



In (3 A 



2 V 2tt 2 

— {(3E c )\ir 1 {f3/i). 
9 



(164) 
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k+1/2 



FIG. 10: Unified scaling diagram of the Coulomb blockade 
in terms of the SET conductance g' and the q . The arrows 
indicate the scaling toward T — (see text). 



The inverse of this quantity is a measure for the width 
Aq of the transition. This width vanishes as T goes to 
zero. 

Notice that Eq. (|164[) is completely analogous to what 
happens at the plateau transitions in the quantum Hall 
regime. In that case we have 



dR H 
dB 



(165) 



i.e. the maximum slope of the Hall conductance Rh with 
varying magnetic field B diverges algebraically in (3. The 
best experimental value of the critical exponent K equals 
0.42. Similar to the experiments on the quantum Hall 
effect, one may consider higher order derivatives of g' 
and q' with respect to <?4£ These quantities diverge even 
faster with higher powers of (3. 

The quantity dq' / dq can directly be measured in the 
experiment since it determines the renormalized gate ca- 
pacitance of the SET C' g = C g dq'/dq. As was shown re- 
cently, the rate of energy dissipation (P) in the SET due 
to a low frequency gate voltage V g {t) = V g + U ul cos cut is 
given by 4 ^ 



fil2 

Poau 2 \Uj 2 9 



9' 



(166) 



Therefore, the maximum in P with varying values of V g 
diverges according to 



r<>2 
9' 



<r 



-(f3E c ) 2 \n-\f3/0. (167) 



Finally, it should be mentioned that the critical behav- 
ior of the SET is likely to change when the effective num- 
ber of channels N^b r ^ between the island and the reser- 
voirs are finite rather than infinite! 14 i 49 Even though we 



expect that our theory of physical obscrvablcs remains 
unchanged, Matveev^ has argued that the critical be- 
havior of the SET can be mapped onto the iV-channcl 
Kondo models This would mean that the transition at 
q' = k + 1/2 becomes a second order one with a finite 
critical value of g' thus closely resembling the more com- 
plicated physics of the quantum Hall effect pi 2 - Progress 
along these lines will be reported elsewhere^ 



B. Quantization of q' 

We have seen that the Thoulcss criterion for the 
Coulomb blockade breaks down at points q = k + 1/2 
where the energy gap A' vanishes. To understand how 
the critical features of the SET permit a change in q' one 
must think in terms of a dynamical process where a unit 
of external charge is added to the system at (imaginary) 
time and removed at r > 0. This process is described 
by the two-point correlation function D(t) given by 



D(t) = ( e »(°)-<*W). 



(168) 



Following Eq. (fTTj) . the tunneling through the SET in- 
volves the sum over all processes D(t) according to 



g' = gT drD(r). 
Jo 



(169) 



From Eq. (fl44j) we obtain the following expression valid 
at T = when q approaches k+1/2 from below 



D(t) = j- 2 e(r)e 



-tA' 



(170) 



This general result includes the isolated island except 
that the AES operators are now renormalized (7 ^ 1) 
and the energy gap A is replaced by the renormalized 
value A'. 

Let us first assume that D(t) denotes the correlation 
of an isolated island. Eq. (| 1 69[) then stands for a semi- 
classical picture of the SET where the island and reser- 
voirs are essentially disconnected. Since the expectation 
value (r) is finite for q < k + 1/2 



(r) = 



f~dTTD(T) 

J °°drD(r) 



1 



(171) 



it is impossible that the tunneling processes described 
by Eq. (|169|) alter the static charge Q on the island. 
However, as one approaches the critical point then the 
expectation (r) diverges. It is thus possible that when 
q passes through k + 1/2, a unit of charge stays behind 
on the island. This extra charge is precisely what lowers 
the energy of the island, i.e. it permits the energy to 
jump from one parabolic branch E c (q — k) 2 to the next 
E c {q-k - l) 2 , see Fig. U 

From the expression for q' in Eq. (|T2"|) it is clear that 
this semi classical picture of the SET gets dramatically 
complicated when the tunneling conductance g is finite. 
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In particular, the second term proportional to g is Eq. 
(f?2")) clearly indicates that the quantization of q' goes 
hand in hand with strong charge fluctuations between 
the island and the reservoirs. Nevertheless, the mecha- 
nism for changing the quasi particle charge q' of the SET 
remains essentially the same. This mechanism solely in- 
volves a vanishing energy gap A'. The only difference 
with the semi classical picture is that the AES operators 
e ±i* m (|168[) generally stand for the quasi-particle 
operators of the SET, rather than those of ordinary elec- 
trons in an isolated island. 

Let us next consider the tunneling process in some 
more detail. We are interested, first of all, in the en- 
ergy difference SE between the states \q + 1} and \q) of 
the SET. Here, \q + 1) is formally defined as follows 



|g+l>= lim \q( T )) 



(172) 



where q(r) = q + 1 for < t < r and q(r) = q for 
tq < r < (3. After elementary algebra we obtain^ 



SE 



A' 



l + e-P A ' 



T 



1 rf 7 2 
d\nT' 



(173) 



Since SE ~ A' at low temperatures (T ^ A') we conclude 
that the transition from \q) to \q + 1) is energetically 
unfavorable. 

Next, the rates for electron tunneling from reservoir to 
island (Tqi) and backward (Tiq) are computed to b o 46 i 48 



<?A' / , A' 
r 0Vio = 4^2 (^coth— ±1 



(174) 



As long as the energy gap A' is finite, the energy dif- 
ference SE in Eq. (|173p and the tunneling rates r i/i 
are not related to one another in any obvious manner. 
However, at the critical point A = we find 



SE = 



1 



r i = Tin = nSE. 



(175) 



Hence, the energy difference between the states \q+l) and 
| q) determines the time the electron resides on the island. 
It is therefore possible that the tunneling processes alter 
the static charge q' of the SET as q passes through fe+1/2. 
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APPENDIX A: ANALYTIC CONTINUATION OF 

K(ioj n ) 

In this appendix we perform the various steps that take 
us from K{iu> n ), Eq. ([50]) . to the expression for K R (u>) 
in Eq. (|84|) . Starting from the correlation function 



we employ the following relations 

de iu)n 



D(iu> m ) 



7T € + lLU m 

de JmD R (e) 



7T e - lLO rn 

and obtain the following expression 

g f deide 2 e 2 lmD R (ei) 



(Al) 

(A2) 
(A3) 



K(iu n ) = 



4~ 



1 



ei + £2 + tu r , 

1 



(A4) 



ei - iuj m e 2 + iuJ m + iuJ n 
Evaluating the sum over u) m we find 

j,,. \ g f de x de 2 e 2 lmD R (e 1 ) 

K(iu n ) = — / — — 

47T J -k z ei - e 2 + iu) n 

x[nb(e 2 ) - n b (ei)]. 

The analytic continuation to real frequencies is now triv- 
ial and we directly obtain Eq. (|84|) . 



(A5) 



APPENDIX B: USEFUL IDENTITIES 

The identities obtained in this Section will be of use 
in Appendix [Cl First, we consider the derivative of the 
average charge Q with respect to g which can be obtained 
as follows 



dQ d 2 lnZ. 



pf 



dg dgdA 
l__d_ 

And A 



Tj2\un\D(iiu n ). (Bl) 



With the help of Eq. (|A3[) . it is convenient to rewrite 
Eq. (|B1[) in the following manner 



where 



UJ n e 



^ — ^ t,i z A- 



(B3) 



A second useful identity for the expression appearing in 
Eq. (|53")) is given by 



Re 



J deI mD R (e) d -^ (B4) 



duj 



9 
'2tt 2 



where the limit u — > is understood. 
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APPENDIX C: EVALUATION OF q TO SECOND 
ORDER IN g 



where the limit u> — > is understood. Next, by expanding 
Eq. (|C5p to the second order in g we finally obtain 



Based on Eqs (|B2j) and (|B4|) we evaluate, in this Ap- 
pendix, the expression for q' in Eq. (|83[) to second order 
in g. We start from the two-point correlation function 
D{iuj n ) which to first order in g is given by 



D(iu> n ) 



tanh(/3A/2) 



9 PY(A) 



7T sinh(/3A) 

^tanh(/3A/2)AfJ^) 
7T u ' 'dA \iu} n - A 

9 \Un\ 



47r (itu n - A) 2 ' 
This result can be written in the following form 
tanh(/3A' 1 /2) 1 g \u 



(CI) 



7i 2 



itu r , 



A' x 4tt (iw„ - A) 2 ' 
(C2) 

where the renormalizcd energy gap and the rcnormaliza- 
tion factor are given as 

A' 1 = A-^(A), t = ^. (C3) 
7r 7^ oA 

Eq. (|C2[) implies the following expression for the retarded 
function: 



D R {e) 



tanh(/3A' 1 /2) 



1 



e-A[+ i0 H 



4?r (e - A + iO+y 



(C4) 



By inserting the result (|C4)) for D R (e) in Eq. (|B4|) we 
find after elementary algebra 



Re 



dui 



g tanh(/3A' 1 /2) dy(Ai) 



2tt 



7i 



a Ax 



Re 



S- tanh 



2tt 

ay 

dA 
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g d 2 y 2 

27 <9A 2 
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(5 BY 2 



9 d ( . d . , 
32^M A M AC ° th 2 



27rsinh/3A <9A 
/3A 



(C6) 



We proceed by inserting the result for D R (e) in the 
expression of Eq. (|B2[) and find 



dQ 
dg 



= 2^9A 

g g 2 
32tt 2 dA 2 



tanh(/3A' 1 /2) 



7? 



y(A;) 



A^coth^) 



(C7) 



(C8) 



Up to second in g the expression for the averaged charge 
Q therefore becomes 



Q = k(q) 



9 



1 

(3Y 2 



D/ 3A 



2tt sinh (3 A 



2tt <9A 
c> 2 



9 



64tt 2 <9A 2 



2 V 2tt dA 



A 2 coth^ ) 



Finally, collecting Eqs. (|C6[) and (|C9[) together we find 
the total result for q' as follows 



9 = Hq) + 



=/3A 



2tt I <9A 



ggy 2 9 9 y 2 9 2 

2tt <9A dA 2tt SA 2 



tanh 



/3A 



_2 ^2 

■^A— ( Acoth 
64tt 2 dA 2 V 2 



V • (cio) 



r 



The result can be written in a slightly more compact 
fashion according to 



Here. 



A' 2 = A - (gMY(A[) 



(C12) 



q' = k{q) + 



1 



1 + eP A 2 



2 q2 
^A—— ( Acoth 

64tt 2 <9A 2 I 2 



9" a d 2 f A ^ 0A 



(Cll) 



represents the second order in g expression for the renor- 
malizcd gap. 

Similarly, with the help of Eqs. (|CT|) and (|C4|) from 
Eqs. (f52"| . (fTTj) . and (|72[) one can compute the other re- 
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sponsc parameters to the second order in g. The results 
can be summarized as follows: 



g'i 



and 

g'u 



in 



g tanh/3Ai/2 
2tT W ' 



(C13) 



k(q) 



1 



g 2 A d 2 
1 + e ^ A 2 ~ 64tt 2 dA 2 



g gAj 

2 7 2 sinh/^A; 
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/3A 



(C14) 



g 

4tt 2 



(C15) 



d f3A 2 ,. / i^A 

TT^ + 0^( Acoth ^ 16) 



/?A 



These results are different from those obtained in the 
leading logarithmic approximation. In some cases (i.e. 



g'jj and q'j) the corrections in g no longer predict an ex- 
ponential dependence on T in the limit where T goes 
to zero. This clearly shows that the expansion to low- 
est orders in g does not provide access to the Coulomb 
blockade phase of the SET. 



APPENDIX D: EVALUATION OF K R {u) 

In this Appendix we present the results of explicit com- 
putations of the response function K R (ui). These results 
can be used, first of all, as an independent check on the 
results of Eqs. fTO5]) . fTTSjl . ijlTSj) . (fHTl) and $TM$>- Sec- 
ondly, they show how the analytic continuation of K(iu n ) 
to real frequencies works in explicit computations. 



1. Weak coupling regime g' ^> 1 



Based on Eqs (|101[) and pi3[) we obtain the following 
expression from Eq. 



J 



iwg 



-in 
2tt 2 



g 2n 2 T g H \ 2itT J 



E 



gE c 

2n 2 



n + 2^T 



e 9/ 2 U^ 



V(i)-V i- 



2ttT J 



cos 2^5 jT^ — 



(Dl) 



Here, the cuttoff n max appears due to the fact that, as usual, we use the low-frequency part of the kernel a(r) only. If 
one takes the proper expression for it into account then one finds n max ~ Ep /T where Ep denotes the Fermi energy. 



2. Strong coupling regime j' < 1 



Given Eq. (|145[) we can write 



„ . , q tanh , 
K R {u>) = - , n n 2 / dee 



n b (e) - n b (A') 
e — A' — uj - i0+ ' 



(D2) 



Hence, 



ImK R (u) 
ReK R (u) 
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47T7 Z 
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47T 2 7^ 



i(A'+w)\ ( iA'\ 2tt /a , n 



By using the following representation of Bosc-Einstcin function in the sum over Matsubara frequencies 

n 6 ( w ) = iw(l + ^)- 1 + 



1 T 

a; 



(D3) 



(D4) 



(D5) 
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wc finally obtain 



K R (w) 
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